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Summary. The Modified Direct Method (MDM) is an iterative mesh smoothing 
method for smoothing planar and surface meshes, which is developed from the 
non-iterative smoothing method originated by Balendran [1]. When smooth planar 
I I meshes, the performance of the MDM is effectively identical to that of Laplacian 

smoothing, for triangular and quadrilateral meshes; however, the MDM outperforms 
Laplacian smoothing for tri-quad meshes. When smooth surface meshes, for trian- 
gular, quadrilateral and quad-dominant mixed meshes, the mean quality(MQ) of all 
mesh elements always increases and the mean square error (MSB) decreases during 
smoothing; For tri-dominant mixed mesh, the quality of triangles always descends 
while that of quads ascends. Test examples show that the MDM is convergent for 
^-H both planar and surface triangular, quadrilateral and tri-quad meshes. 
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1 Introduction 

In finite element analysis it is important always to use high quality meshes: 
low quality meshes lead to unreliable results. A mesh that has been newly 
created usually needs to be improved before it can be used. This improve- 
ment can be made using either (1) mesh clear-up methods, which insert or 
delete nodes as well as change the connectivity of the mesh elements, or (2) 
mesh smoothing methods, which leave the element connectivity unchanged 
and instead reposition the mesh nodes. This paper is mainly concerned with 
smoothing methods. 

There are numerous papers published concerning the topic of mesh smooth- 
ing. In this paper, we just refer some that are the most popular or represen- 
tative of the state of the art for planar and surface meshes. 
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1.1 Related works 

Mesh smoothing methods probably can be classified into four types: the 
geometry-based, the optimization-based, the physics-based and the combined. 
The geometry-based methods obtain new location of nodes by geometric rules, 
local optimization techniques or minimizing objective functions. 

The most popular geometry-based smoothing methods is the Laplacian 
smoothing [12], which repositions each node at the ccntroid of its neighboring 
nodes in one iteration. The popularity of this method comes from its simplicity 
and effectiveness. To improve the basic form of Laplacian smoothing, some 
smart, constrained or weighted variations have been proposed [8,11,27]. 

Another simpler but more effective method is called angle-based approach 
[28], in which new nodal locations are calculated by conforming specific an- 
gle ratios in the surrounding polygons. A geometric element transformation 
method (GETMe) [25] based on a simple geometric transformation can be 
applicable to elements bounded by polygons. 

A projecting/smoothing method is performed for smoothing surface mesh 
by minimizing the mean ratio of all triangles sharing the free node [7]. A 
novel method based on quadric surface fitting, vertex projecting, curvature 
estimating and mesh labeling is applied in biomedical modeling [26] . 

An effective variational method for smoothing surface and volume triangu- 
lations is proposed by Jiao X et al [15], where the discrepancies between actual 
and target elements is reduced by minimizing two energy functions. Also, a 
general-purpose algorithm called the target-matrix paradigm is introduced 
in [17], and can be applied to a wide variety of meshes. 

Different from the geometry-based methods, in optimization-based ones, 
the smoothed position of all nodes is acquired by minimizing a given distortion 
metric. This series of methods is more expensive but can generate better 
results than most geometry-based ones especially at concave regions. Some 
literatures devoted on this topic include [3,10,21,23]. 

The physics-based smoothing methods are the techniques that smooth 
meshes based on physical processing [20], or by solving simple physics prob- 
lems. Shimada [22] proposed a method which treats nodes as the center of 
bubbles and nodal locations are obtained by deforming bubbles with each 
other. A similar algorithm called pliant method is presented in [2]. 

To improve performance, two or more basic methods can be combined into 
a hybrid approach. Some hybrid methods are the combination of Laplacian 
smoothing with optimization-based methods [4,5,9]. 

In the aspect of computation, mesh smoothing methods can be either iter- 
ative or non-iterative. Most of them are iterative algorithms; a non-iterative 
method is proposed by Balendran [1] which is referred to here as the Direct 
Method (DM) . The goal of the DM is simple - to make triangular elements as 
close to equilateral as possible and quadrilateral elements as close to square 
as possible - and it achieves this goal by generating and solving a set of 
optimization equations. 
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1.2 Our contribution 

In this paper we introduce a smoothing method that has the same basic goal 
as the DM, but is iterative rather than non-iterative; we term this method the 
Modified Direct Method (MDM) which can be used to smooth both planar 
and surface meshes. 

The main procedure of MDM for smoothing planar meshes is relatively 
simple: Firstly element stiffness matrices are created based on the type of 
elements. The modified forms of element stiffness matrices are simpler than 
those of DM. And then by assembling all element stiffness matrices, a system 
of Jacobi iteration equations can be formed, which is different from the op- 
timization equations in DM. Finally, the smoothed nodal coordinates can be 
generated by solving the system of Jacobi iteration equations. 

For smoothing surface meshes, the MDM becomes complex: Firstly to 
maintain features of original meshes, feature points are detected and then 
fixed as constrained nodes. Then a Jacobi iteration matrix which is similar to 
that in 2D is also assembled. Thirdly the relocated position of every node in 
each iteration is calculated according to the Jacobi iteration matrix, and then 
projected onto the original mesh to be smoothed position in current iteration. 
And finally, the required smoothed nodal coordinates can be obtained until 
the quality of smoothed mesh will not improve. 

The paper is organized as follows. In Sect. 2, we first give a brief description 
of the basics of the original DM. Then in Sect. 3, we show how the DM 
is modified and developed into the MDM for smoothing planar and surface 
meshes; also we introduce several key techniques such as features detection of 
the MDM for smoothing surface meshes. Finally we present some tests of the 
MDM and make a convergence analysis of the MDM for smoothing surface 
meshes, and summarize them in Sect.4. 

2 The Direct Method (DM) 

The optimization equations solved in the DM are generated by assembling 
element stiffness matrices into a global stiffness matrix. These element stiffness 
matrices arc 6x6 in size for planar triangular elements and 8x8 in size for 
planar quadrilateral elements. The global stiffness matrix for a mesh with n 
nodes is 2n x 2n in size, irrespective of whether the elements are triangles or 
quadrilaterals. 

2.1 Planar triangular mesh 
The element stiffness matrix 

Consider a triangular clement ABC shown in Fig. 1. A* is the position to which 
node A would have to be moved to make the element equilateral, assuming 
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that nodes B and C were fixed; B* is the position to which B would have to 
be moved assuming A and C were fixed; C* is the position to which C would 
have to be moved assuming A and B were fixed. The coordinates of A* are: 



/_ / (1) 



^1 = ^^-^B + Xc) + 2 (i^B + Yc) 




Fig. 1. Triangular element ABC. A*, B*, C* are positions to which nodes would 
have to be moved singly in order to make ABC equilateral. See Eq. 3 



These equations can be rewritten as: 

r 2(Xa - X\) ^ {IXa + OFa) - {Xb + V3Fb) + {-Xc + ^Yc) 

\ 2{Ya - YX) = {QXa + 2Ya) + {V^Xb - Yb) - (VSXc + Yc) ^ ' 

The coordinates of B* and C* are obtained in the same way. The coordinates 
of A, B, C, A*, B* and C* are related by the following equations: 

"2 0-1 -\/3 -1 \/3' 

2 %/3 -1 -Vi -1 
-1^/3 2 -1 -73 

-^3 -10 2 ^/3 -1 
-1 -1^2 

. ^/3 -1 -\/3 -10 2 . 

The left-hand matrix in Eq. 3 is termed the stiffness matrix of the planar 
triangular element. 

The global stiffness matrix and the optimization equations 
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Now assume that ABC is part of a planar triangular mesh that has n nodes. 
Each node of ABC - for instance node A - is then shared with several other 
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elements, and A* can be calculated for each of these. The final position of 
A - its optimal smoothed position - is obtained by averaging the separately 
calculated A*'s (Fig. 2). (This averaging process is effectively identical to that 
used in Laplacian smoothing, which explains why the test results obtained 
using Laplacian smoothing are identical to those obtained using MDM, for 
uniformly triangular and uniformly quadrilateral meshes - see Figs 4, 5.) 
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Fig. 2. Node A belongs to 6 triangular elements. A* can be calculated for each 
element separately (black circles). Optimal smoothed position for A is the centroid 
of these, A* (arrowed open circle) 

The smoothed positions for the complete set of n nodes are given by the 
equations: 
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These equations are the optimization equations that need to be solved in 
the DM. The left-hand matrix - this is termed the global stiffness matrix - is 
created by assembling the individual element stiffness matrices according to 
the connectivity of the elements in the mesh concerned. The elements in the 
global stiffness matrix are obtained during this assembly process, and they 
will of course be different from mesh to mesh. 6^(1 < j < n)is the number of 
elements in the mesh that share node i. 

2.2 Planar quadrilateral mesh 
The element stiffness matrix 

Consider next a quadrilateral element ABCD shown in Fig. 3. A* and C* are 
the positions to which nodes A and C would have to be moved to make the 
element square, assuming B and D were fixed. The coordinates of A* are: 
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Fig. 3. Quadrilateral element ABCD. A*, B*, C*, D* are positions to which diag- 
onally opposite nodes would have to be moved to make ABCD square. See Eq.6 



[ YX = -{-Xb + Xd) + i^{yB+ Yd) 

Similarly as that to triangular element, the coordinates of A, B, C, D, A*, 
B* , C* and D* are related by the following equations: 
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(6) 



The left-hand matrix in Eq. 6 is termed the stiffness matrix of the planar 
quadrilateral element. 



The global stiffness matrix and the optimization equations 

Now assume ABCD is part of a quadrilateral mesh that has n nodes. As before, 
the global stiffness matrix is obtained by assembling the element stiffness 
matrices according to the connectivity of the elements in the mesh concerned. 
The optimization equations associated with this global stiffness matrix are 
identical in form to those for the triangular mesh (Eq. 4). 



3 The Modified Direct Method (MDM) 

The development of the MDM from the DM involves (1) the use of different 
element stiffness matrices, (2) the use of a Jacobian iteration matrix instead of 
a global stiffness matrix, and (3) the replacement of the optimization equations 
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with iteration equations. The mathematical steps involved in this development 
arc broadly similar for the triangular mesh, the quadrilateral mesh and the 
tri-quad mesh. 



3.1 The element stiffness matrices 









1/2 


73/2 








-V^/2 


1/2 


1/2 


-V3/2 








\/3/2 


1/2 








1/2 


V3/2 


1/2 


-73/2 


-73/2 


1/2 


\/3/2 


1/2 









'x\ 




Ya 




Yl 




Xb 




X*B 




Yb 




Y*B 




Xc 




xh 




Yc 







(7) 



Triangular mesh 

According to Eq. 1, the coordinates of A, B, C, A*, B* and C* are then related 
by the following equations: 

1/2 -73/2" 
73/2 1/2 
1/2 73/2 
-73/2 1/2 



These equations can also be rewritten as iteration equations: 

1/2 73/2 1/2 -73/2' 

-73/2 1/2 7^/2 1/2 

1/2 -73/2 1/2 73/2 

73/2 1/2 -73/2 1/2 

1/2 73/2 1/2 -73/2 
-73/2 1/2 73/2 1/2 

X\ and are the coordinates of node A at step A:, and auu 
are the coordinates at step fc+1; the notation for nodes B and C is similar. 
The left-hand matrix in Eq.8 is also the element stiffness matrix for a planar 
triangular element, but it is simpler in form to that used in the DM (Eq.3). 

In Cartesian coordinates, the coordinates of a vertex x. y and z are equiva- 
lent and can be cycled in sequence: X —^Y,Y ^ Z,Z ^ X. Hence, we easily 
extend the stiffness matrix (Eq.8) for a planar triangular element to a surface 
one in 3D (Eq.9). 
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Quadrilateral mesh 



According to Eq.5, The coordinates of A, B, C, D, A* 
also be rewritten as iteration equations: 
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(11) 

The left-hand matrix is also the clement stiffness matrix for a planar 
quadrilateral element, but it is simpler in form to that used in the DM (Eq.6). 
Similarly, the 2D version of stiffness matrix for a planar quadrilateral element 
(Eq.lO) can be extended into 3D version (Eq.ll). 



Tri— quad mesh 

The element stiffness matrices used in a tri-quad mesh are those already given 
for the triangular and quadrilateral meshes (Eq.8, 9, 10, 11). Which of these 
is used for a particular element in a tri-quad mesh depends only on that 
element's type. 

3.2 The MDM for smoothing planar meshes 

Just as in the DM, the element stiffness matrices can be assembled into a 
global matrix, also of size 2n x 2n for a mesh of n nodes. This Jacobi iter- 
ation matrix has several forms, the one shown in Eq.l2 is the most efficient 
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computationally. Therefore it is the one that will need to be solved iteratively, 

starting with the original node coordinates at step and continuing until no 
node needs to be moved by more than the given tolerance distance. 
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, wheyeD = diag{l/ei, 1/ei, l/e2, l/e2, l/e„, l/e„) 

The algorithm for implementing the MDM in 2D has three basic steps: 
(1) Search for elements that share a node, for each node; (2) Assemble of the 
element stiffness matrices into the iteration matrix; (3) Solve of the iteration 
equations until the tolerance distance is reached. 



3.3 The MDM for smoothing surface meshes 

When smooth surface meshes, it is necessary to keep the features of original 
meshes. Many features preservation approaches have been proposed [6,14,16]. 
A popular methods is to classify all the nodes of a mesh into four types: 
boundary, corner, ridges and smooth nodes; and then boundary nodes and 
corner nodes are fixed while smooth node can be adjusted on the whole mesh 
and ridge nodes can only be relocated along the ridges. 

In order to preserve fc^aturcs. we firstly detect the corner nodes and ridge 
nodes via Jiao's approach [14], and fix them as constrained nodes together 
with the boundary nodes although the ridge nodes can be moved along ridges; 
and then after obtaining the relocated positions based on the .Jacobi iteration 
matrix, we project the these new nodes onto the original mesh to get the 
mapped ones which are still the candidates of the resulting nodes in this 
iteration; and thirdly, we check whether there exists inverted elements in the 
incident faces of the mapped node. If does, recover it; otherwise, update the 
target node with the mapped position(Algorithm.l). 

Since the MDM is an iterative smoothing algorithm, we have to conduct 
some indicators to judge when the meshes are smoothed enough and then end 
the iteration. We adopt two mesh quality indicators, the mean quality of mesh 
(MQ) and mean square error (MSE) of all elements, to show target meshes 
are smoothed enough and iterations can break. This can be presented as: 

r MQ'=+^ - MQ'^ < 

otherwise, continue to iterate 

, where Smq and Smse are two user-specified thresholds. Noticeably, the above 
criteria only works when there are no inverted elements. 
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Algorithm 1 The MDM for smoothing surface meshes 

1: Search all iiicklcnt faces F{vi) for each node vt; 

2: Assemble the 3n x 3n Jacobi iteration matrix B = (6ij)o<i<3n,o<3<37i; 
3: while iteration not converge do 
4: Estimate or update normal at each node Vi 
5: for each node Vi do 
6: if Vi is not constrained then 
7: Calculate the relocated node t)f ™ based on B: 

8: Project t;I"''"onto F{vi) to obtain the mapped node 

9: Check inverted element in F(v^'^^). If exists, recover v^""^ •<— Vi 

10: Update all nodes: Vi ^ v^"" 



Normal of vertices 

The direction of each vertex is closely related to its incident faces, for that 
the normal of each vertex is nearly vertical to the normal of any face of its 
incident elements. Therefore, the unit normal of each faces in a mesh should 
be calculated firstly by computing the planar equation of every face. Suppose 
there are m faces shares the vertex u, and the normal of v can be obtained 
by solving the following to x 3 linear equations Nx=l, where N is a to x 3 
matrix whose i*^ row is the unit normal of the i^^ incident face of the vertex 
V, and 1 = (1, 1, is a vector of length to. Since N may be over- or under- 
determined, the' solution is in least squares sense and can be solved by the 
singular value decomposition (SVD). 

Identifying points 

We adopt Jiao's algorithm [14] to detect features based on eigenvalues analysis 
of a symmetric positive semi-definite matrix A: A = N-^WN, where N is a 
TO x 3 matrix we denote in above section, and W be a to x to diagonal matrix 
with Wn equal to the weight. If users do not consider the weights, W can be 
ignored and A is then simplified as N^N. Let Ai, A2 and A3 (Ai > A2 > A3) 
be the three eigenvalues of A. The relative sizes of the eigenvalues Aj of A are 
closely related to the local flatness at a vertex. In general, A has three large 
eigenvalues at a corner, two large ones at a ridge, and one large one at a smooth 
point. Hence, the corners and ridges can be recognized by comparing A3/A1 
and A2/A1 against some thresholds: if A3/A1 > Xo is at corner; if A2/A1 > 
Xr) is on a ridge , where Xc and Xr are two given thresholds. 

Projecting 

After relocating position in an iteration step, smooth node will be projected 
onto its incident faces of the original mesh alone the updated normal. Since 

that, wc just project the; relociatcid point onto its last incident faces (neigh- 
borhood), there maybe not exists a mapped point on the original mesh inside 
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the neighborhood of the above node. Thus, let k be the number of mapped 
points on the incident faces, then wc have: (1) if fc = 0, ignore the relocated 
position and recover its last coordinates; (2) if A: = 1, adopt the only mapped 
point as the new position; (3) if A; > 1, select the nearest mapped point to the 
relocated position as the new position. 



4 Test applications of the MDM 
4.1 Mesh quality 

The simplest way to measure mesh quality is to calculate distortion values for 
each of the mesh elements. The distortion value for a triangular clement should 
measure how close that triangle is to equilateral. One appropriate measure, 
a, was proposed by Lee and Lo [19]. For the triangle ABC shown in Fig.l: 

IICA X CBIl 



= 2a/3- 



CA||V !|AB||V ||BCf 



The value of a lies between and 1; a = when A, B and C are collinear; a 
= 1 when ABC is equilateral. 

The distortion value for a quadrilateral element should measure how close 
that quadrilateral is to square. In this paper we use the measure A proposed 
by Hua [13]; this applies only to convex quadrilaterals. For the quadrilateral 
ABCD shown in Fig.3: 



A 



' ||AB X AD|| • ||BC X BA|| • ||CD x CB|| • ||DA x DC|! 

(||AB||^ + ||AD||2)(||BC||^ + ||BA||^)(||CD||^ + ||CB||^)(||DA||^ + ||DC||^) 



The value of A lies between and 1; A = when any three nodes are collinear, 
i.e., when ABCD is in fact a triangle; A = 1 when ABCD is square. 



4.2 The test applications 

For smoothing planar meshes, a number of test meshes were created and 
smoothed (Figs. 4, 5, 6). The mesh quality results for tri-quad meshes before 
and after smoothing are given in Tables 1. For smoothing surface meshes, a 
tri and a quad surface mesh are created by planar triangulations [18, 19, 24] 
and then interpolated to surface (Fig. 8). 



Table 1. Smoothing for planar tri-quad meshes (T: tri; Q: quad) 



T-dominant 


Q-dominant 


T-MQ T-MSE Q-MQ Q-MSE|T-MQ T-MSE Q-MQ Q-MSE 


Original 0.8688 0.1188 0.8211 0.0713 
LS 0.8756 0.1083 0.8464 0.0571 
MDM 0.8794 0.1083 0.8306 0.0583 


0.9079 0.1243 0.9194 0.0739 
0.9064 0.0936 0.9554 0.0568 
0.9394 0.0737 0.9576 0.0568 
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A tri-dominant mixed mesh is generated by dividing each sliver triangle 
and its neighbors in the original triangular mesh into a smaller triangle and 
a quadrilateral (Fig. 8(e)). Another quad-dominant mixed mesh is created by 
pairing triangles in 2D and then interpolating (Fig.8(g)). 




(a) Original tri mesh (b) Smoothed by LS/MDM 



Fig. 4. Smoothing for planar triangular mesh 

The smoothed surface meshes are shown in Fig. 8. The user-specified 
thresholds £mq and Smse are set as 10~®, 10~^, 10^^ and 10""* for triangu- 
lar, quadrilateral, tri-dominant and quad-dominant meshes, respectively. And 
correspondingly, the MDM converges at the steps 43, 56, 33 and 13. The mesh 
quality results before and after smoothing are given in Table 2. 




(a) Original quad mesh (b) Smoothed by LS/MDM 



Fig. 5. Smoothing for planar quadrilateral mesh 
4.3 Convergence analysis 

In this section, we analyze the convergence of the MDM. When smooth the 
planar meshes, test examples show that the MDM does converge. We have not 
given the mathematical proof for it in theory. For smoothing surface meshes, 
the MDM is much more complicated. As mentioned above, the mean quality 
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of mesh (MQ) and mean square error (MSE) of element qualities are used 
to indicate the quality of whole mesh, wc can calculate and compare the two 
indicators in increasing smoothing iterations to analyze the converge. 




Fig. 6. Smoothing for planar tri-quad meshes 



The triangular surface (Fig.8(a)) is smoothed by MDM in 20, 40, 60, 80 and 
100 iterations. The MQ and MSE of smoothed meshes in different iteration 
stages arc listed in Table 2. And corresponding scatter diagrams are drawn in 
Fig. 7(a) and 7(b). It is clear that the mean quality of mesh ascends while MSE 
descends during increasing iterations. The quadrilateral mesh is smoothed by 
MDM in 10, 20, 30, 40 and 50 iterations. From Table 2 and Fig. 7(c) and 
7(d), we can also receive the same conclusions as that of triangular mesh. 
Similarly, the tri-dominant mixed mesh is smoothed in 10, 20, 30, 40 and 50 
iterations (Table 2, Fig. 7(c) and 7(f)), and the quad-dominant mixed mesh is 
smoothed in 5, 10, 15, 20 and 25 iterations(Table 2). 

According to the scatter diagrams of MQ and MSE, we can learn that: for 
both triangular and quadrilateral surface meshes, the MQ always increases 
and the MSE decreases during smoothing. But the magnitude and rate of 
change are becoming smaller and smaller. Thus, we may in theory receive the 
conclusion that the smoothing iterations will converge after some steps. This 
can be also concluded for quad-dominant mixed mesh. 
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(e) MQs of smoothing tri-domi. mesh (f) MSEs of smoothing tri-domi. mesh 

Fig. 7. Mean quahty(MQ) and mean square error of quaHty(MSE) when smooth 
surface triangular, quadrilateral and tri-domi. meshes. ♦:triangles, •:quadrilterals 



For the tri-dominant mixed mesh, the quahty of triangles always decreases 
while that of quads increases. Noticeably, in beginning iteration steps, both 
of the above qualities change dramatically. Different from the mean qualities, 
the MSE is much more complex: the MSE of triangles significantly ascends in 
beginning steps, and then descends stably; the magnitude and rate of change 
are becoming smaller and smaller. The MSE of quadrilaterals seems to always 
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decline, but obviously, there is an inflection point at steps 30(tlie MSEs at 
both steps 20 and steps 40 arc bigger than that at steps 30); hence, there 
must be a convergence point around steps 30. Our tests (Table 2) proves the 
above conclusion: the MDM converges at steps 33 when the thresholds Smq 
and Emse are set as 10~^. In Table 2, each stage includes 5, 10 or 20 steps. 



Table 2. Qualities for smoothing surface meshes in increasing iteration steps 



Meshes Indicator Stage Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 



Tri mesh 


MQ 


0. 


,8891 


0. 


,8980 


0, 


,9026 


0, 


.9049 


0, 


.9067 


0. 


,9073 


MSB 


0. 


,1565 


0. 


1317 


0, 


,1281 


0, 


.1272 


0, 


.1264 


0. 


,1257 


Quad mesh 


MQ 


0. 


,8592 


0. 


,8814 


0, 


,8929 


0, 


.9027 


0, 


.9122 


0. 


,9198 


MSE 


0. 


,1466 


0. 


,1237 


0. 


.1088 


0. 


.0973 


0, 


.0874 


0. 


,0785 


Tri of tri-domi. 


MQ 


0. 


,8769 


0. 


,8672 


0. 


.8663 


0. 


.8660 


0. 


.8659 


0. 


,8658 


MSE 


0. 


,1509 


0. 


1657 


0, 


.1644 


0, 


.1633 


0, 


.1625 


0. 


,1619 


Quad of tri-domi. 


MQ 
MSE 


0. 
0. 


,7406 
,1397 


0. 
0. 


,7711 
1153 


0, 
0, 


.7750 
.1123 


0, 
0, 


.7772 
.1115 


0, 
0, 


.7785 
.1116 


0. 
0. 


,7792 
,1122 


Tri of quad-domi. 


MQ 
MSE 


0. 
0. 


,9187 
,1066 


0. 
0. 


,9417 
,0884 


0. 
0. 


.9435 
.0839 


0. 
0. 


.9449 
.0813 


0, 
0. 


.9466 
.0777 


0. 
0. 


,9481 
,0747 


Quad of quad-domi. 


MQ 
MSE 


0. 
0. 


,8220 
,1432 


0. 
0. 


,8625 
,1323 


0, 
0. 


.8721 
.1235 


0, 
0. 


.8788 
.1162 


0, 
0, 


.8842 
.1102 


0. 
0. 


,8887 
,1048 



4.4 Tests assessment and summary 

1) For the two topologically uniform planar meshes, i.e., the meshes with the 
same type of element throughout, there is effectively no difference between 
Laplacian smoothing and MDM. (This was referred to earlier, in the descrip- 
tion of the original DM.) Both these methods give smoothed meshes that are 
markedly better than the original unsmoothed meshes. 

2) For the triangular elements in the planar tri-dominant mixed mesh, the 
MDM outperforms Laplacian smoothing while for the quadrilateral elements, 
the reverse is true (Table 1). However, for planar quad-dominant mixed mesh, 
the MDM perfectly outperforms Laplacian smoothing for both triangular and 
quadrilateral elements. 

3) For triangular, quadrilateral and quad-dominant mixed surface meshes, 
the MQ always increases and the MSE decreases in increasing iterations. For 
tri-dominant mixed mesh, the quality of triangles always descends while that 
of quads ascends. 

4) Test examples shows that the MDM is convergent for both planar and 

surface triangular, quadrilateral and tri-quad meshes. 

5) There are some 'sharp' vertices in the smoothed surface meshes. One of 
the probable cause is that both ridge and corner vertices are strictly fixed; the 
other possibly is to project relocated position onto the original meshes rather 
than an uniform underlying surface or a local parametric curve or surface. 
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